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k>" ! Abstract 

IT") ■ Dark Matter (DM) and Dark Energy (DE) can be both described in terms of ultra- 

light Pseudo-Goldstone-Bosons (PGB) with masses tudm ^ 10~ 23 eV and wide ~ 10 _33 eV 
respectively. Following Barbieri et al, we entertain the possibility that a PGB exists with 
mass mi intermediate between these two limits, giving a partial contribution to DM. We 
evaluate the related effects on the power spectrum of the matter density perturbations and 



o 

1-^ \ on the cosmic microwave background and we derive the bounds on the density fraction, // , of 

this intermediate field from current data, with room for a better sensitivity on // in the near 
P-i. future. We also give a simple and unified analytic description of the free streaming effects 

<-j I both for an ultra-light scalar and for a massive neutrino. 

> : 

k*n : 1 Introduction and statement of the problem 

To understand the nature of Dark Matter (DM) and Dark Energy (DE) is one of the most intriguing 
problems in all of today's physics. Proposals for possible solutions are not lacking. More difficult 
is to prove any of them by a direct signal; an example in the case of DM would be the detection 
in the laboratory of a WIMP interaction. In this note we pursue the idea that part of DM may be 
due to an ultra-light Pseudo-Goldstone-Boson (PGB), arising from the spontaneous breaking at a 
scale close to M Pi of an extended approximate symmetry, and we study a possible related signal. 

The idea that DM or DE may be interpreted as the energy density of a PGB is not new. 
The axion is the prototype example for DM. Well known is also the fact that the mass of the 
hypothetical scalar field associated to DE would have to satisfy the bound m DE ^ 10 _33 ey. The 
tightness of this bound provides in fact a clear case for interpreting also this scalar as a PGB [U El Ej. 
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Following this view, a calculable microscopic model for the potential of a PGB for DE has 
been proposed in Ref. |1] , based on an approximate U{l) n flavour symmetry of the right-handed- 
neutrinos. For the purposes of the present paper, the key feature of the model is the presence of 
several PGBs, Gi, up to five in the general version, with a potential of the form 

V(Gi) = E if 4(co8(G i /F i ) + l), (1) 

where F { are mass scales related to the scales of spontaneous breaking of the U(l) factors and 
the heights of the various potential terms, are controlled by small parameters that explicitly 
break the U{l) n flavour symmetry. For the interpretation of DE, one such term, with i = DE, 
must have G EE « F DE ~ M Pi and ^de ~ 3 ■ 10~ 3 eV, with G l £ E the initial value of the DE- 
PGB field. On the other hand, the presence of other similar terms in eq. (0), with i ^ DE, 
makes it natural to ask what the manifestation of the other PGBs could be. If their masses, 
rrii ?s yU^/Fj, are greater than m DE « H , the Hubble constant today, then these extra PGBs 
contribute to DM, with possible effects on the growth of cosmic clustering. One such scalar, of 
mass m E M ~ 1CT 23 -j- 10~ 22 eV, has in fact already been invoked with the purpose of suppressing 
the apparently unobserved cusps arising in the halos of conventional Cold Dark Matter (CDM) 
0. In the potential of eq. (0), the parameters that would lead to the interpretation of DM in 
terms of such a field are G l EM ~ F E m ~ O.OlMpi and Udm ~ 

30 eV. 

Here we entertain the possibility that a PGB exists with a mass m I intermediate between 
™dm ~ 10~ 23 eV^ and tude ~ 10~ 33 eV. In the range of parameters 



M 



pi 



0.01^1; & = 3 • 1(T 3 -r- 30eV (2) 



one would have indeed 



mj = 1(T 33 + lO~ 23 eV; // = ^- ~ 10~ 3 -r- 1 (3) 



where // is the fraction of the energy density in the field Gj relative to the matter density, possibly 
given by a heavier PGB. In general, for G i p ~ Fj, it is, when fj < 1, 

fi ^i^?™<^W-) l/2 ) (4) 
6 Mpi n eq 

The interest of the presence of such a field is that it would lead to peculiar features of the 
matter power spectrum, quite analogous to but also generally different from the ones of a massive 
neutrino. As we now recall, these effects are related to the wave properties of the fluid associated 
with an ultra-light scalar, which inhibit the clustering of its energy density at sufficiently small 
scales. The detection of these effects might therefore constitute a signal of the overall underlying 
picture. 

The CMB is also affected by such a field. In fact, when the expansion rate H is higher than 
mj, the field freezes and behaves as a cosmological constant. If this occurs at the scale factor a F 
before the present, the field behaves for a < a F as a cosmological constant with density fraction 
f2/<(a) ~ (a/a P ) 3 fiQ m (assuming matter dominates). If f2/<(a) is not negligible at decoupling 
time, then there is a distortion in the anisotropy spectrum due to the early Sachs- Wolfe effect. 
Since decoupling occurs at Hd w 6 ■ 10~ 29 eV, we expect the CMB distortion will be maximal for 
m; near this value. 
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2 Free streaming in an analytic approximation 



The physics underlying these phenomena is well known EH- A straightforward modification 
of the CMBFAST software [B] indeed allows a numerical calculation of these effects (See Sect. 3). 
To ease the reading and to make explicit their connection with the effects of a massive neutrino, 
we derive them in a simple and effective analytic approximation. 

What makes the difference between an ultra-light scalar or a massive neutrino and ordi- 
nary CDM is the presence of a pressure term in the Euler equation for the equivalent fluids, 
(/c 2 /a 2 )Cg^5fc, where 8 k is the density perturbation of the mode with comoving wave number k, a 
is the scale factor and c e // is an effective sound speed, which becomes less then 1 at sufficiently 
large a. For a PGB of mass it is 

C ft Rj for a > , (5) 

n 2a m<t> ~ 2m/ v ; 

whereas for a neutrino of mass m u ^ HH] 

for a > — — = clnr (6) 



m,,a ~ m,, 



where T° is the neutrino temperature today. This pressure term introduces a Jeans length, 1/kj, 
below which the fluids are unable to cluster: they "free stream". Equating the pressure term to 
the source term, 4nGp8k, during matter domination, one finds 

k*(a) = LSea^XWn^ and fcfta) = 1.22a^H Q^. (7) 

A classic result is that, if all matter contributing to the cosmic density is able to cluster, like 
ordinary CDM, the corresponding density fluctuations, 8%, all grow like the scale factor, 8% oc a, 
between matter-radiation equality, a = a eq , and now, a = a . Similarly, if only a fraction (1 — /) 
can cluster, the growth is slower |TT] 



5 c k <xa p ,p= 1/4(^/1 + 24(1 (8) 

This means that, relative to the case in which free streaming is neglected, the power spectrum of 
the matter density fluctuations is given by 

r(k) = Ptrue = ( Qmaa: )2(p-i) (9) 

Pno— free— streaming ^min 

where (a m j n , a max ) is the /c-dependent interval of scale factor during which free streaming is rele- 
vant. 

In the case of the PGB, defining a m as the scale factor at which the field starts to oscillate, 
at H ~ m<£, and dj(k) as the scale factor at which k = kj(a), it is a min = max(a m ,a eq ) and 
o-max = min(a ,aj(k)). The relevant scales are therefore (m 30 = m^lO -30 ^)) 

kf q = k%a eq ) « 9 • W^mlfMpc- 1 (10) 
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(oo) « 1.1 ■ lO-^C^^^Mpc- 1 (11) 
fc£ = fcj(a m ) « 3.6 • lO-W^ 3 (h 2 a m Y^Mpc- 1 (12) 



and, r'P, unaffected below k^ in , for k > fc^ in becomes 



r*(Jfc) « (7f^) 8(1-p) (13) 



where fc^ in = max(kf qJ fc*) and fc£ aa; = mm(fc, fctf). 

In the neutrino case, one proceeds in a fully analogous way. Let us restrict ourselves to the 
case m v < leV, so that the neutrino becomes non relativistic after matter-radiation equilibrium, 
o>nr> a>eq- In this case there are only two relevant scales (m eV = m u /eV) 

k v NR = k»(a NR ) « 0.026m e 1 (, 2 (/i 2 a m ) 1 / 2 Mpc- 1 (14) 

kl = k v j(a ) « 2.2m ey (/i 2 fi m ) 1/2 Mp C - 1 (15) 



and, for > /catr, 



r"(*) - (^)^ (16) 



mux 



where k v max = min(k, k^) 

Eq.s (fT^ITH|l are the reference formulae. They are immediately useful when the power spectrum 
with neglect of free streaming is simple to compute, since only then, from eq. (JSJ), Ptrue is easily 
determined from r(k). This is so if the fraction /, relative to conventional CDM, of the extra 
component we are considering is small, below 10%, which is anyhow the case of interest here. 
If so, r^(k) in eq. (fl~3|) and r u {k) in eq. (IT6| give directly the power spectrum normalized to 
conventional ACDM with 3 massless neutrinos. In principle, this is not precise. In the low mass 
range, at sufficiently early times but still after equilibrium, neither the PGB nor the neutrinos 
behave as CDM, even apart from free streaming. Since massless neutrinos are anyhow included, 
this is a negligible effect for neutrinos. In the PGB case, a better approximation of the denominator 
in eq. (0), especially at small scales, consists indeed in including the contribution of a fictitious 
scalar field with c e // = 0. 

Note that in the neutrino case, its mass determines as well /, which appears in the exponent 
of eq. ([HI), according to 

^ L m 'v iL m 

Note also that, if more than one neutrino is massive, all below 1 eV so that each relative / is 
small, one obtains the power spectrum normalized to the massless case by taking a factor as in 
eq. (|T6|) for every massive neutrino. 

Figs. 1 and 2 show the power spectrum of density perturbations for the PGB (Fig. 1) and the 
massive neutrino cases (Fig. 2), normalized as described above, for some representative values of 
the relevant parameters. The continuous lines are from the numerical code, whereas the dashed 
lines represent the analytic expressions discussed above. The analogies and the differences between 
the PGB and the neutrino case are apparent from the figures themselves and from the discussion 
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above. Neutrino masses are known to exist at a level that might be within reach of future 
cosmological measurements. In the ideal case in which all the various parameters of ACDM 
were known, there are two characteristic differences between Fig. 1 and 2, that can be used to 
distinguish them. The interval in k with a non- vanishing slope is controlled by the mass parameters 
in a characteristically different way for the neutrino and the PGB cases. The slope itself depends 
on the fraction /, which in the PGB case is a free parameter whereas for the neutrinos it is tied 
to the mass by eq. ifTTj) . 

3 Comparison with observations 

We now compare with current observations the matter power spectrum and CMB when the extra 
PGB scalar is present along with a cosmological constant and ordinary CDM. Both the cosmologi- 
cal constant and ordinary CDM could be effective descriptions of ultra-light PGBs, of appropriate 
mass, but this is irrelevant to the discussion of this Section. As recalled in the previous Section, 
the perturbation equations for a generic scalar field can be expressed as a perfect fluid with a 
scale-dependent sound speed. Moreover, we take the field behaving as a cosmological constant 
(wi = —1) for H > mi and as dark matter (wi = 0) afterward. 

We produce the power spectra and the temperature and polarization CMB spectra by using a 
suitably modified version of CMBFAST [E] which includes the extra fluid. We explore the range 
of masses and densities 

mj = (1(T 6 -r 10 8 )m 30 
// = 0-1 

which includes as extreme cases the behavior of a cosmological constant (masses smaller than 
H « lCT 33 eV ) and of ordinary cold dark matter (masses larger than the typical galaxy size, 
10" 24 eV). 

We compare the results via a grid-based likelihood analysis with four observations: the power 
spectrum of SDSS [12] ( cutting at k — 0.2Mpc//i); the power spectrum of Lyman-a clouds [T3| ; 
the variance cr 8 = 0.90 ± 0.03 derived by [TH] combining Lyman-a data in the SDSS catalog with 
WMAP data. Finally, we use the CMB temperature and polarization spectra of the 2001 WMAP 
dataset [H]. Notice that the datasets are all independent except for the use of WMAP to fix 
the initial perturbation amplitude in deriving <t 8 : when we combine all constraints we neglect 
this residual correlation. (It is also to be noted that as is subject to considerable systematic 
uncertainty [T5] . potentially quite larger than the statistical one.) For the SDSS we used the 
likelihood routine provided by M. Tegmark which includes the full covariance matrix while for 
CMB data we adopt the likelihood routine provided by the WMAP team [TJ]. When comparing to 
SDSS and Lyman-a spectra we marginalize over the amplitudes (separately for the two datasets), 
so that the information on the spectrum amplitude is used only in the third test. In this way the 
large-scale structure constraints are well separated into those derived from the large scale slope 
(SDSS), the small scale slope (Lyman-a) and the amplitude (cr 8 ). Finally, we combine all four 
tests to give an exclusion plot for mj, fj. 
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We let the other cosmological parameters vary in the ranges 

G 0.6^0.8 
G 0.06^0.15 
G 0^0.3 
G 0.95 -=-1.05 

while to save computer time we fix the baryon fraction f2&/i 2 = 0.023. We assume flat priors for 
all parameters. 

The results are shown in Figs. 3 (the likelihood for each dataset) and 4 (the combined likeli- 
hood), after marginalization over the other cosmological parameters (and the overall normalization 
for WMAP, SDSS and Lyman-a). The Lyman-a constraints turn out to be the strongest. As 
expected, there is an intermediate-mass window of observability between 10~ 31 eV and 10~ 23 eV. In 
this region of interest for the mass mj, the relative density // is clearly limited to be below 10%, 
which is a significant constraint but leaves ample room for a relevant component of a DM-PGB, 
as given in eq. ((3J). 

The prospects for detecting the signal from an intermediate-mass PGB scalar are interesting. 
The same cosmological observations that will constrain the neutrino mass can provide limits to the 
PGB mass and abundance, since the PGB scalar affects both the growth of fluctuations and the 
power spectrum shape. In particular, deep weak lensing observations are expected to provide in 
few years upper limits to the neutrino abundance f v of less than one percent ^Hl, mostly through 
the damping effect on the fluctuation growth. Since a scalar with the same abundance // provides 
a very similar fluctuation growth (at least for mj G 10~ 30 -j- 10 _24 eV) we can expect that also // 
can be constrained to better than the percent accuracy. This similarity of behavior raises in fact 
also the issue of the level of degeneracy of scalars with neutrinos, a problem which may deserve 
attention. 
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Figure 1: Comparison of the ratios r v {k) of numerical matter power spectra (red full curves) with 
the analytic approximation discussed in the text (green dashed curves). Upper panel: 3 massive 
neutrinos with mass 1/3 eV each, and £7a = 0,0.3,0.7, (ordered top to bottom according to the 
left plateau). Lower panel: same, with masses 1/9 eV. 
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Figure 2: Comparison of the ratios r^(k) of numerical matter power spectra (red full curves) 
with the analytic approximation discussed in the text (green dashed curves). Upper panel: PGB 
with fraction // = 0.05 and masses m 30 = 0.1, 1, 10 (ordered top to bottom according to the left 
plateau), and £7a = 0. Lower panel: same, with = 0.7. 
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Figure 3: Likelihood functions at 68,95 and 99.7% c.l. (dark to light gray) for the parameters 
m 30 = m I /10~ 30 eV and fj obtained marginalizing over r, Q m h 2 , n s and h, and fixing VL b h 2 = 0.023. 
The date employed are discussed in the text. The discreteness of the grid causes some wiggling 
in the function. 
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